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Abstract 


The constrained Dirichlet boundary value problem x = f{t,x), 2;(0) = 
x(T), is studied in billiard spaces, where impacts occur in boundary 
points. Therefore we develop the research on impulsive Dirichlet prob¬ 
lems with state-dependent impulses. Inspiring simple examples lead to 
an approach enabling to obtain both the existence and multiplicity re¬ 
sults in one dimensional billiards. Several observations concerning the 
multidimensional case are also given. 
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1 Introduction 

For three hundred and fifty years several important problems of mechanics 
and physics have been investigated and modeled as mathematical problems with 
impacts. Starting with the works of G. D. Birkhoff, dynamical systems in spaces 
of the billiard type have been intensively studied. The simplest impact law, for 
absolutely elastic impacts, can be described geometrically as the equality of 
angles before and after a collision with a boundary of the billiard space. This 
law, for simplicity, will be assumed in the present paper. 

An elementary observation is that the dynamical system of a billiard type 
with a uniform motion can be modeled by the simple impulsive second-order 
system 



( 1 ) 


where K = intK C M" is a compact subset, and I is an impulse function 
describing the impact law. It is easy to check that for a unit ball B{0, 1) C M" 
and for the equality of the angle of incidence and angle of reflection, one has 
I{x{s),x{s)) = —2(x{s),x{s))x{s). 

Let us imagine a one-dimensional billiard which is not a straight line but 
a graph of some differential function 7 : [a, b] —>■ K. We can think about 
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some hills and valleys on our simple one-dimensional table. Assume the grav¬ 
ity directed downstairs. Then, the horizontal component of the acceleration 
is nonzero. In fact, the motion can be described by a more general equation 
x{t) = —fcgrad 7 (x(t)), where k is some constant depending on the gravity. 

If we allow an external force depending on time (e.g., a wind), we get even 
more general equation x{t) = f{t,x{t)). One can also easily guess that for 
tables generated by a nondifferentiable but Lipschitz function 7 we obtain a 
second-order differential inclusion x{t) € F{t,x{t)). This all above motivates to 
investigate the system 


x(t) = f{t,x{t)), 

x(s-l-) = x(s) -I- /(x(s), x(s)), 


for a.e. t > 0,x{t) G intK 
if x(s) G dK. 


( 2 ) 


Different kinds of problems in billiards are interesting, e.g., existence of 
periodic motions and their stability, number of distinct periodic trajectories, 
etc (see [ 6 ] and references therein). We are interested in the boundary value 
problem of the Dirichlet type: [@ and x(0) = x(T) = 0] motivated by the 
research explained below. 

In last decades a lot of papers on impulsive boundary value problems have 
been published. Most of them concerns impulses at fixed moments. In this case 
one can see direct and clear analogies with the approach and results for prob¬ 
lems without impulses. Several difficulties appear when impulses depend on the 
state variable or both the time and state. One meets the case in e.g. differential 
population models or mechanics where impulses occurs if some quantities attain 
a suitable barrier. The papers dealing with state-dependent impulsive problems 
focus attention mainly on initial or periodic problems, and results on the exis¬ 
tence, asymptotic behavior or a stability of solutions. We refer e.g. to [T] (and 
references therein), where the impulsive periodic problem, also for some kinds 
of second-order differential equations, is studied. 

Unfortunately, the Dirichlet impulsive boundary value problem cannot be 
brought to the first-order one, and different techniques are needed (comp, some 
recent papers using a variational approach for problems with fixed impulse times, 
[?, ?, ?]. Quite recently in [S] the authors examined the problem 

xit) = f(t,x(t)), for a.e. t G [0,r], 

x(s-l-) = x(s) -I- /(x(s)), ifs = g(x(s)), ^ 

where g is a. function satisfying some additional conditions. They provided 
a new method to solve the problem. Namely, they successfully transformed the 
problem to the fixed point problem in an appropriate function space. Note that 
both the impulse function I and the barrier Graph{g) = {(x, s);s = g{x)} are 
not adequate for billiard problems. Indeed, in billiards the impulse depends 
also on the velocity before the impact. Moreover, the barrier is not a graph of 
a function with arguments in a phase space. Note also that in [5] the assump¬ 
tions insist all trajectories go through the barrier without coming back while in 
billiards trajectories stay on the same side of the barrier after the impact. The 
technique presented in [ 8 ] does not work for the Dirichlet problem in a billiard 
space. This problem, as far as the author knows, is still unexplored. 
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Therefore the aim of the paper is to study the existence and multiplicity of 
solutions to the Dirichlet impulsive boundary value problem 

( x{t) = f{t,x(t)), for a.e. t S [0,T],x{t) G intK, 

< i;(s+) = i;(s) + J(a;(s), i(s)), ii x{s) G dK, (4) 

[ a^(O) = x{T) = 0, 

where / describes the impact law of 

(HO) the equality of the angle of incidence and angle of reflection and the equal¬ 
ity of a length of the velocity vector before and after the impact. 

By a solution of (jH) we mean a continuous function x : [0, T] —>■ R", which has an 
absolutely continuous derivative in intervals where x{t) G intK, satisfies x{t) = 
f{t,x{t)) for a.e. t G [0,r], the impulse condition i;(s-|-) = i:(s) -|- I{x{s),x{s)) 
in boundary points, and the boundary value condition x(0) = x(T) = 0. Note 
that different solutions of problem (I5|) can have different numbers of impact 
points. In general they also can slide along the barrier dK, and their integral 
representation is difficult. So, the problem attracts by its nontriviality. 

We start in Section [5] with two examples inspiring the further research. Then 
the results about a one-dimensional case are presented in Section [3] The final 
section concerns problems in R" so in multidimensional billiard spaces. The 
author hopes the paper is a valuable development of the area of second-order 
impulsive boundary value problems. Simultaneously, the theory of billiards is 
developed in the direction of crooked tables. 

2 Two inspiring examples 

Let K = [—a, a], where a > 0, and let G be the Green function for the 
autonomous problem x = 0,a:(0) = x{T) = 0, see e.g. [9]. We assume in the 
whole paper that 

(HI) the right-hand side / (comp. (|3])) is a Caratheodory function, shortly 
/ G C'ar([0,T] x R), which is, for simplicity, integrably bounded, i.e., 
fi-,x) : [0,T] —>• R is measurable for every x G R, /(t, •) : R —?> R is 
continuous for a.e. t G [0,r], and |/(t,x)| < m{t) for some integrable 
function m G L^([0,T]) and each (t,x) G [0,T] x R. 

Then solutions of the nonimpulsive Dirichlet problem can be obtained as fixed 
points of the operator T in (^^([OjT]), 

J'(x)(t) := / G{t,s)f{s,x{s))ds. 

Jo 

Assume that x G FixT. Then |x(t)| < M for some M > 0 and every t G [0,T]. 
If M < a, then the solution x remains in intK and the problem (|3]) is trivially 
solved. If M > a, then the trajectory bounces off the barrier, and the situation 
becomes nontrivial. 

In the first example below we consider simple autonomous equation. 


3 


Example 2.1. Consider the problem 


r x{t) = 2, for a.e. t G [0,1], |a;(t)| < i, 

< i;(s+) = -x{s), if \x{s)\ = i, (5) 

[ x(0) = a:(l) = 0. 

It means that I{x{s),x{s)) = —2x{s). 

Notice that the unique solution of the corresponding nonimpulsive problem 
takes the form x{t) = t{t — 1) with ||a::|| = 1/4. Therefore, it does not fit in 
K = [—1/8,1/8]. The first time when the trajectory meets the barrier can 
be easily computed as tp = ■ After the impact, we have a trajectory 

x[t) = + {^/2 — \)t + ■ One can check that, after some impacts, we 

obtain x{l) ^ 0. To find a solution of (O we notice, at first, that the unique 
solution of the autonomous nonimpulsive problem 

/ x{t) = f{x{t)), forte [0,1], , . 

\ x{Q) = cc(l) = 0 

satisfies the equality a;(l/2—t) = x(l/2+t) for every t e [0,1/2]. Indeed, one can 
check that y(t) := a;(l — t) is a solution, and, by the uniqueness, x(t) = a;(l — t). 

Now, we would like to find an initial velocity v = x(0) such that the function 
starting as x{t) := t{t+v) has the following moments of impacts: {ti, 1/2,1 —ti} 
and x(l) = 0, so we look for a solution of ([5|) symmetric with respect to t = 1/2. 
It is sufficient to find v such that 4ti + 2to = 1, where tp is a negative solution of 
the equation t{t + z;) = 1/8 (see Figure[T|). After some computation we obtain 
V « -0.8568, to « -0.127 and h « 0.186475. 



> 


t 


Figure 1. Solution of ([5]) for x(0) ~ —0.8568. 

The solution of problem is not unique. To find the second one it is 
sufficient to solve the equation 8ti — 6tp = 1 (we obtain the initial velocity 
vi ~ —1.76579). It has 7 impact points. In fact, there are infinitely many 
solutions obtained when the absolute value of the initial velocity tends to infinity. 

The second example concerns a nonautonomous case, so we loose a symmetry 
with respect to t = T/2. 
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Example 2.2. Consider the problem 

f xit) = 6t, for a.e. t £ [0,11, bit)! < 

\ i;(s+) = -x{s), if b(s)| = |, (7) 

[ x(0) = a:(l) = 0. 

Now, the function x{t) = — t is the unique solution of the corresponding 

nonimpulsive problem. But the barrier is such that x(l/2) = —3/8 and i:(l/2) < 
0. Hence, we have an impact. In fact, it is the only impact before t = 1. Denote 
by xi the function which equals x up to the impact, xi{t) = 6t for t >ti = 1/2, 
and lims^i /2 ^i(s) = —d:i(l/2). One can check that xi{t) = — j for 

t > 1/2. Therefore, a::i(l) = 1/4 (see Figure[2]). 



Figure 2. Solution of x{t) = 6t, a;(0) = x{l) = 0, and the trajectory with the impulse effect. 

The initial velocity of xi is equal to a;i(0) = —1. We know that each solution x^, 
depending on the initial velocity x«(0) = u, of the impulsive Cauchy problem 

f xit) = 6t, for a.e. t G [0,11, |x(t)| < |, 

J x(s+) = -x(s), if |x(s)| = |, (8) 

[ x(0) = 0 

up to the first impact has the form Xy{t) = + vt. It is not hard to check 

that for every u G [—1,1] we get x„(l) > 0, and for — -^243/256 < u < 1 the 
trajectory meets only the upper part of the barrier (x = 3/8) (see Figure [3|). 
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Taking into account the above trajectories, we try to solve problem © by 
increasing the absolute value of v. 

Assume that u < —1. The first impact time ti is a solution of = —vti—3/8. 
After this impact the trajectory is given by x{t) = C — {dtl +u)t — 4wti —9/4 (we 
use an impact law in ti). Analogously, after the second impact time t 2 we obtain 
x{t) = — { 6 t 2 — 6t1 — v)t + 3/8 + 5 t 2 — Gt‘lt 2 — vt 2 - Now we evaluate in t = 1 and 

ask for v such that x(l) < 0. Since the last point a;(l) depends continuously 
on V (see details in Section [3]) , it follows that there exists a trajectory with 
a;(l) = 0, so a solution of ©. One can check that for v = —1.218 we have 
a;(l) Ri -0.0006379. 

The above two examples show that we can expect solutions of © for big 
initial velocities, and we can also expect a multiplicity of solutions. 


3 Existence and multiplicity of solutions 

We start with a simple consequence of assumption (HI). 

Proposition 3.1. Assume that x : [0,T] —> AT = [—a, a] is a solution of the 
Cauchy problem 

x{f) = f{t,x{t)), for a.e. t > 0, 

x{s+) = x{s) + I{x{s),x{s)), ifx{s)GdK, (9) 

a;(0) = 0, 

with an initial velocity x(0) = v, and • ■ • the moments of impacts. Then, 
for each sufficiently big |?;|, one has x{tk)x(tk+i) < 0 for every k >1, that is, 
each two consecutive impact points are on the opposite components of the barrier 
dK. 

Proof. It is sufficient to check that |i(f)| > d for some d > 0 and every t G [0, T], 
if V is sufficiently big. But we know, that 

x{s)ds = V + / f{s,x{s))ds, 

Jo 
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up to the first impact time ti > 0, so 

\x{t)\>\v\- f \f{s,x{s))\ds>\v\- f m{s)ds. 

Jo Jo 

Furthermore, |i(ti+)| = |i(ti)|, and 

x{t) = ±{ti+) + f f{s,x{s))ds 
Jti 

for t G (ti,t 2 ], where t 2 is the second impact time. Hence, for t G [ 0 ,^ 2 ], 

|i(t)| > |i;(ti+)| — f m{s)ds > — f m{s)ds. 

•J t\ 0 

We proceed up to T, and obtain 

\x{t)\ >d:= |?;| - ||m||i > 0, 

for every t G [0,T] and |u| > ||m||i, where ||m||i is an L^-norm of m. □ 

To simplify our considerations, in the rest of the section we will assume that 
(H2) the right-hand side / is Lipschitz with respect to the second variable, i.e., 
\f{t, x) — f{t, y)\ < a{t)\x — y\ for a.e. t G [0, T] and some a G T^([0, T]). 

We are in a position to formulate the following 

Theorem 3.2. Let f satisfy (HI) and (H2), and let Xy,Xw be two solutions 
of ® with initial velocities v,w > ||to||i, and impulse times ti,...,tk and 
si,..., Sk+j, respectively. If j > 2, then there exists a solution of (|4]). 

The key point in the proof of this theorem can be formulated as 

Lemma 3.3. Let f satisfy (HI) and (H2), and letb > ||m||i. Then the operator 
Vt : R —t [—r,r], Vt{v) := Xy{T), is continuous in b. 

Proof. Assume that Xb has k impulse times ti(&),... ,tk{b), and tk{b) T, i.e., 
Xb{T) G (—r,r). The case |a:f,(T)| = r, very similar, is considered at the end of 
the proof. In the sequel, by X(^s,p,v) we will denote the unique solution of the 
nonimpulsive Cauchy problem 

f x{t) = f{t,x{t)), forte [0,T], 

I x{s)=p, (10) 

[ x(s) = V. 

Step 1. We show that 

for every p > Q there exists <5 > 0 such that, for h' G {b — 6,b + d), 

the solution Xb' has exactly k impacts, and \ti{b') — ti(b)\ < p for any (11) 

t e {l,...,fc}. 
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Take any 77 > 0 so small that {ti{b) — rj, ti{b) + r])r\ {ti+i{b) — rj, U+i{b)+r]) = 0 
for i = 1,..., fc — 1, and assume, without any loss of generality, that Xb{tk{b)) = 
—r. Then (see Proposition IXTD Xb{t) < 0 for every t G {tk-i{b),tk{b)), and 
Xb{tk-i{b)) = r. Since Xb{tk-i{b)+) < 0, one has X(tfc_i(6),r,±(tfc_i(b)+))(0 < 
—r for any t > tk{b) close to tk{b). Take t G {tk{h),tk{b) + ij) such that 
^(tfc_i(b),r,i(tfc_i(f))+))(f) < —X- By the continuous dependence on initial con¬ 
ditions for (ITOl) . there exists 0 < 5i(fc) < 1 such that the inequalities 


\s - tk-i{b)\ < 6 i{k), \p-r\<Si{k), \v - Xb{tk-i{b)+)\ < Si{k) (12) 


imply that the solution X(^s,p,v) satisfies X(^s,p,v) (i) < Hence, it has an impact 
in some tk < t < tk{b) + r]. 

On the other hand, Xb{[tk-i{b)^tk{h) — ry]) C (—r,r] which implies that 
there exists 52 {k) > 0 such that 52 {k) < 6 i{k) and inequalities (IT^ imply that 
I ~ X(^tk-i{b),r,xi,{tk-i(b)+)) I |c'-([0,T]) < V S-Ild X(^s,p,v) ([^fe-1 (^), tk{b) — f]]) G 

(-r, 00 ). Thus tk > tk{h) - p. 

Now we start to analyze the behavior of solutions near the point (tk-i(b), r). 
Choose 0 < (53(fc) < 1 such that 


h{k) < 


S 2 {k) 


ftk-i{b)+S3{k) 


2(||a;h|| -I-1) 


and 2 


Jtk-i(b)-S3{k) 


m{s)ds < 


S 2 {k) 


(13) 


Analogously as above we find 62 {k — 1) > 0 such that the inequalities 


\s-tk-2{b)\ < S 2 {k- 1 ), \p+r\ < S 2 {k- 1 ), \v-Xb{tk-2{b)+)\ < S 2 {k- 1 ) (14) 


imply that ^(tfe_ 2 (b),—r,i 5 (tfc_ 2 (b)-t-)) I Ic^([o,t]) ^ ^sik') and has an 

impact in some G itk-iib) — 63 {k),tk-i{b)+ 63 {k)). Then, after the impact, 
we obtain a function denoted by x such that x(t) < 0 and, since d 3 {k) < 1, 

|i(T)| < |ib('r)| + 1 < Ilibll -I- 1, 


for every r G {tk-i,tk-i{b) + 63 {k)). 
Consequently, 


r > x{tk-i{b) + Ssik)) = x{tk-i)+ 

ptk-i{b)+S3(k) 


ftk-i{b)+S3(k) 


x{T)dT 


ifc-1 


> r — 


|i(r)|dT > r - 253(fc)(||i:t,|| -|- 1) > r - 52 {k). 
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Moreover, assuming without any loss of generality that tk-i < tk-i{b), 


\x{tk-i{b) + Ssik)) - Xb{tk-iib) + ^ 3 (fc))l = 

|•tk-l{b)+S3(k) 

= x{tk-i+)+ / /(r,i(r))dT 

Jtk-i 

^tfc_i(b)+i53(fe) 

-Xb{tk-l{b) + )- f{T,Xb{T))dT 

ptk-iib)+S3{k) 

-x{tk-i)+ /(r,i(r))dT 

ptk-i{b)-\-S3{k) 

+Xb{tk-l{b))- f{T,Xb{T))dT 

J tfc_i(b) 


l•tk-l{b)+S3{k) 


' tfc-i 

rtk-i{b) 


f{T,x{T))dT 

|•tk-l{b)+S3(k) 
k-i(b) 


-x{tk-l) + 

r^k-iK^) r^k- 

+Xb{tk-l)+ / f{T,Xb{T))dT - / f{T,Xb{T))dT 

'^ik-1 '^ik-1 

ptk-lib) 

< \x{tk-i) - Xb{tk-i)\ + / |/(T,i(r)) +/(r,Xb(T))|dr 

ptk-iib)-\-S3{k) 

+ \fiT,x{T)) - f{T,Xb{T))\dT 

Jtk-i{b) 

|•tk-l{b)+S 3 {k) /\n(k) dn(k] 

< \x{tk-i) - Xb{tk-i)\ +2 m{T)dT < H-— = 52{k). 

Jtk-i ^ ^ 


We proceed, and find, at last, 5 := <52(1) G (0,1) such that ||a;(o_o,b')~2^(o_o,6)llci([o,T]) < 
< 53 ( 2 ) for any 6 ' S (5 — <52(1), b + <52(1)), and X(o,o,b') has the first impact ti in 
(ti(b) - < 53 ( 2 ),ti(6) + < 53 ( 2 )), where 


Ssik) < 


S2(k) 

2(l|ib|| + 1) 


and 


ftk-i(b)+S3(k) 


I tk-i{b)-S 3 {k) 


m{s)ds < 


<52 (fc) 
2 


Then, for b' G {b — 62 ( 1 ), b + <52(1)), the solution Xb' of problem (0 has exactly 
k impacts, and \ti{b') — ti{b)\ < 77 for any i G {1,... ,k}. 

A proof for Xb{tk{b)) = r is the same. 

Step 2. Take an arbitrary e > 0. Analogously as in Step 1 we assume that 
Xb{tk{b)) = —r, and we choose 0 < 0 < 1 such that the inequalities 


\s-tkib)\<9, |p + r|<6», \v - Xb{tkib)+)\ < 9 


imply that - X(^tk{b)-r,xbitk(b)+))\\ < £■ In consequence, |(T) - 

X{tic(b), — r,xi,(tkib)+)){T)\\ < £■ 
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Choose 0 < T] < 9 such that 


Q rtk{0)+r) 0 

^ < oni ■ II , 1 ^ and 2 / m{s)ds < 

and apply m to find (5 > 0. Now, like in Step 1, we can check that, for every 
b' G {b — 6, b + 6), the solution Xf satisfies 

\tkib') - tkib)\ < 9, \xb'(tk{b')+) - Xbitkib)+)\ < 9. 

Thus \xb'{T) — Xb{T)\ < e. 

To finish the proof of the lemma it is sufficient to notice that the case tk (b) = 
T is quite similar. The only difference is that, for b' close to b, the solution Xb> 
can have only k — 1 impulses. But still this solution considered on a little larger 
interval, e.g. [0,T + 1], must have an impulse near T, and the above proof 
arguments work. □ 

Remark 3.4. One can prove the continuity of Vr without the assumption 
b > ||m,||i, because of the continuity of the impulse function I{z) = —2z, but a 
formulation of Lemma 13.31 is sufficient for our considerations. 

Now we are able to prove Theorem 13.21 

Proof. We will show that there exists u{X) = (1 — X)v + Xw (A G [0,1]) such 
that the solution Xu{\) of ([9]) with the initial velocity u(A) solves (jl]). 

Without any loss of generality, assume that Xv{tk) = r. Take A ^ 1. Then 
Xu{\){T) attains the barrier twice (on opposite sides), at least, because new 
impacts occur on different parts of the barrier (see Proposition 13.11) . From the 
Darboux theorem, since Xu{\){T) depends continuously on A (see Lemma 13.311 . 
there exists A G [0,1] such that Xu{x) {T) = 0. The proof is complete. □ 

We finish our consideration with the observation that, indeed, there are 
initial velocities v and w satisfying assumptions of the above theorem. 

Theorem 3.5. Let f satisfy (HI) and (H2), v ^ 0, and Xy be a solution of ([3]) 
with k impulse times. Then there exists c > 0 such that the solution Xcv of 
has at least k + 2 impulse times. 

Proof. Let ti{v),..., tk{v) be impulse times for Xy. Divide the interval [0, ti(r’)] 
into three intervals of equal length. 

At first, we denote by Xy and Xcv the solutions of a nonimpulsive Cauchy 
problem x = f{t, x),x{0) = 0, and show that there is cq > 0 such that \Xcvit2) — 
Xcv{ti)\ > 2,r for c > cq and for each < 1,^2 G [0,T] with \t 2 — ti| > ti[v)/'i. 
Indeed, if cv is an initial velocity, then 

lira(i)I > lira(0)1 - / |m(s)|(is > c|z;| - ||m||i, 
io 
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for every t G [0,r], see the proof of Proposition 13.11 Take ci > 0 such that 
ci\v\ — ||m||i > 0, and let c > ci. Then 

1 

Wxcvih) - Xcv{ti)\ = / \xcv{s)\ds> {c\v\-\\m\\i)-ti{v). 

Jti J 

It is sufficient to find cq > ci > 0 such that §(co|n| — ||TO||i)ti(z;) > 2r. It is 
easy to see that cq > 

Now, if c > Co, then the solution Xcv has an impulse time ti{cv) in [0, ti{v)/3). 
Furthermore, we have the second impulse time t 2 {cv) in [ti(cv),ti(cv)+ti{v)/?>) 
and the third one t^[cv) G \t 2 {cv),t 2 [cv) +ti{v)/3). Thus t^{cv) < ti{v). 

We can choose cq so big that 

\xcvit)\ > \xv{t)\ for every t G [0,T]. (15) 

Indeed, we know that \xcv(t)\ > c|z;| — ||m||i for every t G [0,T]. Similarly we 
check that \xv{t)\ < |f| + ||”^||i for every t G [0,T]. So, it is sufficient to have 
co|u| — ||m||i > |u| + ||m||i which gives cq > (2||m||i + |c|)/|u|. 

Without any loss of generality we can assume that Xcv (t) > Xv {t) > 0 up to 
the first impulse time ti{cv). Obviously, 

rtiiv) rtiiv) 

Xcv{ti{v)) = / Xcv{s)ds> / Xy{s)ds = r, 

Jo Jo 

hence, Xcv(ti(cv)) = r (the upper part x = a oi the barrier). From the previous 
considerations we know that Xcv(t 2 {cv)) = —r and Xcv{io{cv)) = r. From (1151) 
it follows that 

t^+i{cv) < ti+i{v) for every j > 0. 

Indeed, for i = 1, if ^(cc) < ti(c), then ti(cv) < t 2 {v). If ti^cv) > ti(u), then 
Xcv(ti(v)) G (—r,r), and Xcv(ti(v)+) and Xy(ti(v)+) are negative. Moreover, 
Xcvit) < Xvit) for t > ti{v) up to the next impulse time. Thus, again, ti(cv) < 
t2{v). For i > 1 we proceed analogously. In consequence, ti+ 2 (cw) < ti(v). 

□ 

Besides giving the existence of solutions of (IH, the above theorem implies 
their multiplicity, as we can see in the following. 

Theorem 3.6. If f satisfies (HI) and (H2), then there exist infinitely many 
solutions of O- 

Proof. One can easily modify a proof of Theorem 13.51 to obtain a solution of 
problem (|9]) with 2k impulse times before ti{v) for each fc > 1. For each k we 
get at least k different solutions of problem dU, and the proof is complete. 

□ 
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4 Multidimensional billiards. Perspectives of re¬ 
search 


A multidimensional case is much more complicated, and only some special 
situations are investigated below. We give ideas how one could solve the problem 
and indicated some difficulties. Some ideas are formulated as open problems. 
Hence, we treat the section as an impulse to do deeper research. We start, like 
in the previous section, with an inspiring example. 

Example 4.1. Let the billiard space be equal K = B{0,r) C the closed 
ball with a center 0 and radius r > 0. Consider the Dirichlet problem 

{ x{t) = a := ( 01 , 02 ), for a.e. t G [0,T], 

i(s+) = a;(s) +/(a;(s), i:(s)), if x(s) G dK, (16) 

cr(0) = x(T) = 0, 

where I{x{s),x{s)) = —^{x{s),x{s))x{s). Therefore the line 02 X — aiy = 0 is 
invariant under the flow generated by x{t) = a,x(0) = 0,i;(0) = (aoi,ao 2 ). It 
implies that we can solve the problem (HU) as one-dimensional (see Example 
12.11) . Hence, there are infinitely many solutions of (ITH)) . 

The above example leads to the following result for simple billiards (with a 
uniform motion). 

Theorem 4.2. Assume that K C K", n = 2k, is a compact smooth (by ‘smooth’ 
we mean, here and in the sequel, at least of class ) manifold with boundary, 
0 G intK, and K is strongly star-shaped with respect to 0, i.e., tx G intK for 
every x £ K and t G [0,1) (see, e.g., p. 11 for a definition of star-shaped 
sets). Then the Dirichlet problem 

( x{t) = 0, for a.e. t £ [0,T], 

< x{s-\-) = x{s)-\-I{x{.s),x{s)), ifx{s)£dK, (17) 

[ a;(0) = x{T) = 0, 

with I satisfying the standard impact low (HO) (see Section\^ has a nontrivial 
solution. 


Proof. It is sufficient to find a ray {tz; t > 0}, for some z G OK, such that the 
vector Oz is perpendicular to the tangent hyperspace to dK at z. Indeed, if z 
satisfies this, then the function x{t) = vt with v := ^z attains dK at z and in 
the time T/2, then changes the velocity to 


X 



1 -^ 





= -V, 


and finally arrives at 0 = x(T). 

To prove the existence of z mentioned above, we take a small ball B{0, r) C 
intK and denote by R : dB{f), r) —dK a projection along rays, existing by the 
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strong star-shape of the set K. For each point y G dK we take an outer normal 
vector n{y) to K at y and denote by N(y) its projection onto the tangent 
hyperspace to 5'"“^ at y/\y\. We define the map g : 5'"“^ —>■ K", g{w) := 
N{R{rw)). Thus we have obtained a tangent vector field on the sphere 5"“^. 
From the hairy ball theorem (see [3], Cor. 7.4, p. 238) it follows that g{w) = 0, 
for some w G which means that z := R{rw) is the vector we required.. □ 

Let us remark that a regularity of K may be slightly weakened. For instance, 
if dK is a manifold, then the map JV(-) is Lipschitz and, moreover, the set K 
is an absolute neighborhood retract (see, e.g., [3], p. 500). The latter property 
will be used in a proof of Theorem 14.81 

Before the second example let us describe the situation we will deal with. 
Let K C M" be a compact smooth manifold with boundary, and 0 G intK. 
Consider problem (|4]) with I describing the standard impact law (HO) and / 
satisfying 

(Hl)„ / : [0, r] X K” —>• R” is an integrably bounded Caratheodory function, i.e., 
f{-,x) : [0,T] —^ R” is measurable for every x G R”, /(t, •) : R" —>■ R" 
is continuous for a.e. t G [0,r], and \f{t,x)\ < m{t) for some integrable 
function m G L^([0,T]) and each {t,x) G [0,T] x R", 

and (H2). Notice that / can be given explicitly as I{y,v) := —2{v,n{y))n{y), 
where, for some s > 0, y = x{s) G dK, v = x(s) and n{y) is an outer normal 
vector to dK at y. Obviously, I is continuous. Smoothness of the set dK 
together with the continuity of I justifies formulating the following 

Hypothesis 4.3. The operator Vr ■ R" —t K, Vt{v) := Xy{T), where Xy is a 
solution of the problem , is continuous. 

On the way to prove it one meets some difficulties. Some solutions of problem 
([2|) can slide along the boundary dK, at least for some time, and then their 
dynamics is essentially changed, i.e., x{t) ^ f{t,x{t)) on some set with nonzero 
measure. This can happen if {x{s),n{x{s))) = 0 in a boundary point a;(s) G dK. 
Therefore we have a nontrivial hybrid system, where the dynamics on dK can 
be described not only by / but also by a regularity of dK (by a normal vector 
to dK). We hope one could describe it explicitly if Lf is a sublevel set of some 
smooth function. We leave the proof of Hypothesis 14.31 to further studies. 

In what follows we show how one can use the above hypothesis. 

Example 4.4. Assume that K, f and I are as above and define, for every 
d > 0, an attainable set in T as Ad := {Xy(T)-, |n| = d}. By Hypothesis 14.31 it is 
a hypersurface (not necessarily smooth) in K. Let be the area surrounded 
by Ad. Assume that 0 G fid and 0 ^ flad for some a > 0 (see Figure |4]). Then, 
from the continuity of Vr it follows that there exists c = (1 — X)a -I- A, A G [0,1], 
such that 0 G Acd. This means that problem (|3]) has a solution with an initial 
velocity |no| = cd. 
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Figure 4. Attainable sets for juj equal to d,ad,cd. 

The above example leads us to the idea of using the winding number of a 
curve, in a 2-dimensional case, or, more generally, a topological degree in R" 
(see, e.g., [7], p. 6-8, where the degree is presented as the rotation of the vector 
field). 

Indeed, as above we define the sets = {{a;„(T); \v\ = d} C K C K” 

and Qd C K, the area surrounded by A^. Assume that 0 ^ A^. Obviously, 
All can be treated as the image A^ = where —>■ R” is 

defined as V^{v) := Xdv{T). Since is continuous, the topological degree 
deg(V"y, 5(0,1)) is well defined. Moreover, the family {V,^}d>o is continuous, 
by Hypothesis 14.31 Thus we get to 

Theorem 4.5. Let f and I satisfy (Hl)n-(H2) and (HO), respectively. Assume 
the truth of Hypothesis \4.3\ and that 0 ^ Ad for some d > 0. Then 

(i) deg(V^,5(0,1)) ^ 0 implies the existence of a solution Xy to problem 

with an initial velocity v € B{0,d), 

(ii) deg(V^b 1)) deg(V^^, 5(0,1)), for some di < d 2 , implies the 
existence of a solution Xy to problem m with an initial velocity v G 
5(0,d2)\5(0,di). 

Proof. The first part easily follows from the existence property of the degree. 
The second one from the homotopy invariance. Indeed, if there were no zero of 
Vfi for each di < d < d 2 , then, since 14^^ and are homotopic, one would 
have deg(Hyb 5(0,1)) = deg(Hyb 5(0,1)); a contradiction. □ 

A natural hypothesis is that fast movements, i.e., big velocities, can compen¬ 
sate the unevenness of the billiard table described by the nontrivial acceleration 
f{t,x). To check this, we note the following 
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Lemma 4.6. Under the assumptions (Hl)n-(H2) and (HO), for every e > 0 
there exists d > 0 such that, for each |z)| > d, each solution Xy of (ED, with an 
initial velocity v, and for each t G [0, ti(v)] (ti{v) is the first time of impact, see 
Section\^ one has \x(t) — vt\ < e. In other words, the movement is close to the 
one with a uniform velocity. 

Proof. Take a small p > 0 such that p||m||i < min{e, 1}. Put D := 5{K) + 1, 
where 5{K) is the diameter of the compact set K. Now, let d > 0 be such that 
dp > D. This implies that the curve z{t) = vt, for |t;| > d, attains the boundary 
dK before the time p, and dist{z{p),dK) > 1. Then, for a solution X(o_o,„) of 
the Cauchy problem x = f{t, x), x(0) = 0, i(0) = v we obtain 

|a:(o,o.'»)(0 - 2:(i)| < / |i:(o,o.«)(s) - < tllmlli < £ 

Jo 

for every t < p. Moreover, |a;(o,o,«)(p) ~ ^{p)\ ^ Pll”^lli < 1 which implies that 
^(o.o.«)(p) ^ H and, consequently, a;(o,o,j,) meets the barrier dK at a time before 
p. This completes the proof. □ 

The problem is that the increase of |u| causes the increase of a number 
of impacts. Nevertheless, since the impulse function / is continuous, we can 
formulate the following hypothesis which is left as an open problem. 

Hypothesis 4.7. Under the assumptions (HI)n-(H2) and (HO), for every e > 0 
there exists d > 0 such that for each |u| > d, each solution Xy of ED, with an 
initial velocity v, and for each t G [0,T], one has \x{t) — vt\ < e. 

Now we are in a position to formulate a result on a correspondence between 
the general problem Q and the uniform motion Dirichlet problem (ITTl) . i.e., the 
Dirichlet problem in the Birkhoff billiard space. 

Theorem 4.8. We hold the assumptions of Theorem |^.5[ assume the truth 
of Hypothesis and assume that there is a sequence dk oo such that 
deg(W^'“, i3(0,1)) 0 and deg{W^’‘~^,B{0,1)) ^ deg(W^'“, 5(0,1)), where 

: S‘^~^ —>■ M”, := Zy{T) and Zy is a unique solution of the problem 

HID with the initial velocity v. Assume also that dist(W^*'(S'" ^),0) > tq > 0 
for some ro and every fc > 1. 

Then there exists ko > 1 such that deg{V^*‘, B{0, 1)) ^ 0 for every k > ko. 
In consequence, problem O has infinitely many solutions. 

Proof. Since K is an ANR (see, e.g., [2]), take £ > 0 such that each £-near maps 
from S”“^ to K are ro-homotopic (see [5], p. 111). Choose d > 0 for £ as in 
Proposition 14.71 and ko such that dkg > d. Then dist(V'y''(S”“^), 0) > 0 for 
every k > ko, and, by the homotopy property of the degree, 

deg(Py^^ 5(0,1)) = deg(<^5(0,l)) y^O. 

Moreover, by Theorem 14.51 (ii), there is a solution Xk with an initial velocity 
lufel G (dfe_i,dfc). □ 
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Example 4.9. Take K := 5(0, r) C M", where n is odd, and T := 1. Then 
= ids^-^ and = —zdsn-i. It implies that deg(tTy^“^,5(0,1)) = 

1 and deg(IT^^“^, 5(0,1)) = —1. Moreover, dist(IE^^“^(5"“^), 0) = 1. Thus, 
all geometric assumptions from Theorem 14.81 are satisfied, and the theorem 
works for any integrably bounded right-hand side /. 

Open problems and perspectives of research 

1 The author realizes that a computation of deg(Ey,5(0,1)) is difficult in 
general. He leaves as an open problem the question on sufficient conditions 
for deg(H^,5(0,1)) ^ 0. It will be a subject of a further study. Let us 
note that it is not hard to check that, if deg(IEy,5(0,1)) ^ 0 for some 
d > 0, then for small ||m||i one has deg(Hy, 5(0,1)) ^ 0, as well. Hence 
problems with a uniform motion should be carefully examined. 

2 Is it true that limsup^_^oo I deg(Hy, 5(0,1))| > 0 for billiards which are 
strongly star-shaped with respect to 0? Is it true for strictly convex bil¬ 
liards? 

3 For simplicity we have assumed that the right-hand side / implies the 
uniqueness of a solution of problem ® . It is interesting and quite natural 
to consider multivalued right-hand sides because of the control theory, 
where we ask for a suitable control u under which the Dirichlet problem 




a:(s-|-) = x{s) -I- I{x{s), i(s)), if x(s) G dK, 
a;(0) = x{T) = 0, 


has a solution. Indeed, we can define F{t,x) := {f{t,x,u)-,u G U{t)}, 
where U(t) is an admissible set of controls, and examine the inclusion 
X G F{t, x). 

4 It is worth checking the results presented in the paper under weaker growth 
conditions on the right-hand side / (comp, assumption (HI)„). 

5 In the whole paper we have dealt with impulse functions satisfying condi¬ 
tion (HO), that is, with fully elastic collisions with barrier. The non fully 
elastic case is left for the future study. 
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